The Central Idea

A chemical reaction model usually combines three kinds of structure:

A compact modeling architecture is

\[ \boxed{ \text{composition} +\text{reaction mechanism} +\text{rate law} +\text{transport} +\text{energy} \longrightarrow \text{system evolution} }. \]

Chemical Species as State Variables

A reaction model tracks quantities associated with chemical species. Common state variables include

The best choice depends on whether the system is closed or open, homogeneous or spatially distributed, and whether volume or pressure changes materially.

Stoichiometry

Consider the general reaction

\[ \boxed{ aA+bB\rightarrow cC+dD }. \]

The coefficients \(a,b,c,d\) determine the proportional amounts in which reactants are consumed and products are formed.

Stoichiometry constrains quantities even before a kinetic model is specified.

Stoichiometric Coefficients

It is often convenient to assign signed stoichiometric coefficients

\[ \nu_A=-a,\qquad \nu_B=-b,\qquad \nu_C=c,\qquad \nu_D=d. \]

Reactants receive negative coefficients and products positive coefficients.

This sign convention makes reaction-network equations compact.

Extent of Reaction

The extent of reaction \(\xi\) provides one way to represent progress of a single reaction:

\[ \boxed{ n_i=n_{i,0}+\nu_i\xi }. \]

Every species changes according to the same reaction progress variable scaled by its stoichiometric coefficient.

Conservation of Atoms

Chemical reactions rearrange atoms rather than creating or destroying them. A valid reaction model must therefore conserve each chemical element unless the model explicitly includes nuclear processes.

Element balances provide a powerful verification test for stoichiometric models.

Mass Conservation

For an ordinary closed chemical reaction system,

\[ \boxed{\frac{dm_{\text{total}}}{dt}=0}. \]

In an open system, total mass changes only through inflow and outflow:

\[ \boxed{ \frac{dm}{dt} =\dot m_{\text{in}}-\dot m_{\text{out}} }. \]

Reaction changes composition, but not total mass.

Concentration

Molar concentration is

\[ \boxed{C_i=\frac{n_i}{V}}. \]

If volume is constant, concentration changes directly with mole number. If volume varies, dilution or compression must also be included.

Reaction Rate

Reaction rate measures the speed of chemical transformation. For

\[ aA+bB\rightarrow cC+dD, \]

a consistent rate \(r\) can be defined through

\[ \boxed{ r =-\frac{1}{a}\frac{dC_A}{dt} =-\frac{1}{b}\frac{dC_B}{dt} =\frac{1}{c}\frac{dC_C}{dt} =\frac{1}{d}\frac{dC_D}{dt} } \]

when the chosen concentration formulation and volume assumptions make these relationships appropriate.

Rate Laws

A rate law relates reaction rate to state variables and parameters. A general empirical form is

\[ \boxed{ r=kC_A^mC_B^n }. \]

The exponents \(m\) and \(n\) are reaction orders. They must be determined from mechanism or data; they are not automatically equal to the stoichiometric coefficients of the overall reaction.

Mass-Action Kinetics

For an elementary reaction,

\[ aA+bB\rightarrow\text{products}, \]

mass-action kinetics gives

\[ \boxed{ r=kC_A^aC_B^b }. \]

This relationship follows from the assumed elementary reaction mechanism, not merely from balancing the overall chemical equation.

Zero-Order Reactions

A zero-order rate is independent of reactant concentration:

\[ \boxed{ -\frac{dC_A}{dt}=k }. \]

Integrating gives

\[ \boxed{ C_A(t)=C_{A,0}-kt }. \]

First-Order Reactions

A first-order decay model is

\[ \boxed{ -\frac{dC_A}{dt}=kC_A }. \]

Its solution is

\[ \boxed{ C_A(t)=C_{A,0}e^{-kt} }. \]

This exponential structure appears in radioactive decay, simple decomposition, pharmacokinetic elimination, and many effective reaction models.

Half-Life of a First-Order Reaction

Setting \(C_A=C_{A,0}/2\) gives

\[ \boxed{ t_{1/2}=\frac{\ln 2}{k} }. \]

The first-order half-life is independent of initial concentration.

Second-Order Reactions

A simple second-order model is

\[ \boxed{ -\frac{dC_A}{dt}=kC_A^2 }. \]

Integrating gives

\[ \boxed{ \frac{1}{C_A} =\frac{1}{C_{A,0}}+kt }. \]

Unlike first-order kinetics, the half-life depends on initial concentration.

Temperature Dependence

Reaction rates commonly increase with temperature. The Arrhenius equation is

\[ \boxed{ k=Ae^{-E_a/(RT)} }. \]

Here \(A\) is a pre-exponential factor, \(E_a\) is activation energy, \(R\) is the gas constant, and \(T\) is absolute temperature.

Linearized Arrhenius Form

Taking logarithms gives

\[ \boxed{ \ln k =\ln A-\frac{E_a}{R}\frac{1}{T} }. \]

A plot of \(\ln k\) against \(1/T\) is linear under the Arrhenius model, with slope

\[ \boxed{-\frac{E_a}{R}}. \]

Reversible Reactions

A reversible reaction

\[ A\rightleftharpoons B \]

can be modeled with forward and reverse rates:

\[ r_f=k_fC_A, \qquad r_r=k_rC_B. \]

The net rate is

\[ \boxed{ r_{\text{net}} =k_fC_A-k_rC_B }. \]

Chemical Equilibrium

At equilibrium, forward and reverse processes balance so that macroscopic composition stops changing:

\[ \boxed{ r_f=r_r }. \]

Equilibrium does not mean molecular reactions stop. It means opposing processes balance.

Equilibrium Constants

For the idealized reaction

\[ aA+bB\rightleftharpoons cC+dD, \]

a concentration-based equilibrium expression is

\[ \boxed{ K_c= \frac{C_C^cC_D^d} {C_A^aC_B^b} }. \]

More rigorous thermodynamic models use activities rather than raw concentrations.

Reaction Quotient

The reaction quotient \(Q\) has the same algebraic form as the equilibrium expression but is evaluated using the current state.

Comparing \(Q\) with \(K\) indicates which direction of reaction is thermodynamically favored under the assumptions of the model.

Kinetics and Thermodynamics Are Different

Thermodynamics addresses which states are favored and what equilibrium is possible. Kinetics addresses how quickly the system moves.

A reaction can be thermodynamically favorable yet proceed extremely slowly because of a kinetic barrier.

Consecutive Reactions

A reaction sequence

\[ A\xrightarrow{k_1}B\xrightarrow{k_2}C \]

gives the coupled system

\[ \boxed{ \frac{dC_A}{dt}=-k_1C_A } \] \[ \boxed{ \frac{dC_B}{dt} =k_1C_A-k_2C_B } \] \[ \boxed{ \frac{dC_C}{dt}=k_2C_B }. \]

Intermediate \(B\) can rise and then fall because it is both produced and consumed.

Parallel Reactions

If

\[ A\xrightarrow{k_1}B \qquad\text{and}\qquad A\xrightarrow{k_2}C, \]

then

\[ \boxed{ \frac{dC_A}{dt} =-(k_1+k_2)C_A }. \]

Competing pathways determine selectivity between products.

Reaction Networks

Large chemical systems contain many species and reactions. Let \(\mathbf{C}\) be the vector of species concentrations and let \(\mathbf{r}(\mathbf{C},T)\) contain reaction rates.

A compact homogeneous-system model is

\[ \boxed{ \frac{d\mathbf{C}}{dt} =S\mathbf{r}(\mathbf{C},T) }. \]

The matrix \(S\) is the stoichiometric matrix.

The Stoichiometric Matrix

Each column of \(S\) corresponds to a reaction and each row to a species. Entries record the signed stoichiometric change in that species when the reaction advances.

Matrix form separates network structure from kinetic formulas and is valuable for large reaction systems.

Conservation Laws from Stoichiometry

Linear combinations of species may remain invariant under the reaction network. These invariants correspond to conservation relationships and can often be identified from the stoichiometric matrix.

They provide both conceptual insight and numerical verification checks.

Reaction Mechanisms

An observed overall reaction may result from several elementary steps. A mechanism proposes those steps and the intermediates connecting them.

Mechanistic modeling can explain rate laws that cannot be inferred directly from the balanced overall equation.

Steady-State Approximation

If an intermediate \(I\) is produced and consumed rapidly, its concentration may remain approximately constant over a relevant interval:

\[ \boxed{ \frac{dC_I}{dt}\approx0 }. \]

This approximation can reduce a complex mechanism to a simpler effective rate law.

Rate-Limiting Steps

A reaction network may contain one or more slow processes that control the overall time scale.

Identifying a rate-limiting step can support model reduction, but it should be justified quantitatively rather than assumed from intuition alone.

Catalysis

A catalyst changes reaction pathways and kinetic barriers without being consumed in the net reaction.

Models of catalysis may include adsorption, desorption, surface coverage, catalytic-site occupancy, and competing reaction pathways.

Enzyme Kinetics

A common effective model for an enzyme-catalyzed reaction is the Michaelis-Menten equation:

\[ \boxed{ v= \frac{V_{\max}S} {K_M+S} }. \]

Here \(S\) is substrate concentration, \(V_{\max}\) is the limiting rate, and \(K_M\) sets the concentration scale over which saturation becomes important.

Saturation

When \(S\ll K_M\),

\[ v\approx\frac{V_{\max}}{K_M}S, \]

so the response is approximately linear.

When \(S\gg K_M\),

\[ \boxed{v\approx V_{\max}}. \]

The system becomes saturated.

Batch Reactors

A batch reactor contains a fixed amount of reacting material during the reaction period, with no continuous inflow or outflow.

For constant volume,

\[ \boxed{ \frac{dC_i}{dt} =\sum_j \nu_{ij}r_j }. \]

Batch models are direct applications of reaction-network ordinary differential equations.

Continuous Stirred-Tank Reactors

An ideal continuous stirred-tank reactor, or CSTR, assumes perfect mixing. For species \(A\),

\[ \boxed{ V\frac{dC_A}{dt} =F C_{A,\text{in}} -F C_A +V R_A }. \]

Here \(F\) is volumetric flow rate and \(R_A\) is the net rate of production of \(A\) per unit volume.

Residence Time

For a reactor of volume \(V\) with volumetric flow \(F\),

\[ \boxed{ \tau=\frac{V}{F} }. \]

Residence time is a characteristic time scale for continuous-flow systems.

Plug-Flow Reactors

A plug-flow reactor allows composition to vary along the reactor direction. Under suitable steady assumptions, a species balance may take the form

\[ \boxed{ F\frac{dC_A}{dx} =A_cR_A }. \]

The state evolves with position rather than with time for the steady model.

Conversion

Conversion measures the fraction of a reactant consumed:

\[ \boxed{ X_A= \frac{n_{A,0}-n_A}{n_{A,0}} }. \]

Conversion is dimensionless and often more useful for reactor design than raw concentration alone.

Yield and Selectivity

When multiple products are possible, engineers may care about how much desired product is produced relative to reactant consumed or undesired product formed.

Definitions of yield and selectivity vary by context, so they must be stated explicitly in the model.

Energy Balances

Reaction and temperature are often coupled. A simplified reactor energy balance may be written

\[ \boxed{ \rho C_pV\frac{dT}{dt} = \dot Q +(-\Delta H_r)V r +\text{flow terms} }. \]

Exothermic reactions release heat, while endothermic reactions absorb heat.

Thermal Feedback

Temperature can increase the rate constant through the Arrhenius relationship, while reaction changes temperature through heat release or absorption:

\[ T\longrightarrow k(T) \longrightarrow r \longrightarrow \dot Q_{\text{reaction}} \longrightarrow T. \]

This feedback can create multiple steady states, ignition, extinction, or thermal runaway in some reactor models.

Multiple Steady States

Nonlinear chemical reactor models may possess more than one steady solution for the same external operating conditions.

Identifying all physically meaningful steady states and their stability is crucial when operating near transitions.

Stability of Chemical Systems

Suppose

\[ \frac{d\mathbf{x}}{dt} =\mathbf{f}(\mathbf{x}). \]

Near an equilibrium \(\mathbf{x}^*\), local stability can be studied through the Jacobian matrix

\[ \boxed{ J= \left[ \frac{\partial f_i}{\partial x_j} \right]_{\mathbf{x}=\mathbf{x}^*} }. \]

Eigenvalues of \(J\) determine local growth or decay of small perturbations.

Oscillating Chemical Reactions

Some nonlinear reaction systems exhibit oscillations rather than monotonic relaxation toward equilibrium.

Such systems demonstrate that chemical dynamics can contain feedback, instability, periodic behavior, and pattern formation.

Diffusion

Molecular transport by diffusion is commonly modeled with Fick's law:

\[ \boxed{ \mathbf{J}_i=-D_i\nabla C_i }. \]

The negative sign indicates movement down a concentration gradient in the simplest model.

Diffusion Equation

Conservation combined with Fick's law gives

\[ \boxed{ \frac{\partial C}{\partial t} =D\nabla^2C }. \]

This model describes spatial smoothing of concentration in a homogeneous medium.

Reaction-Diffusion Systems

Chemical reaction and diffusion combine as

\[ \boxed{ \frac{\partial C_i}{\partial t} = D_i\nabla^2C_i +R_i(\mathbf{C},T) }. \]

These models appear in combustion, catalysis, electrochemistry, atmospheric chemistry, biological pattern formation, and materials processing.

Advection-Reaction-Diffusion

When a flowing fluid transports species,

\[ \boxed{ \frac{\partial C_i}{\partial t} +\mathbf{v}\cdot\nabla C_i = D_i\nabla^2C_i +R_i }. \]

This couples chemistry to fluid mechanics and transport.

Characteristic Chemical Time

For a first-order process,

\[ \boxed{ \tau_{\text{chem}}=\frac{1}{k} }. \]

Comparing chemical time with transport or residence time helps determine whether reaction or transport controls system behavior.

The Damköhler Number

A Damköhler number compares a transport or flow time scale with a chemical reaction time scale:

\[ \boxed{ \mathrm{Da} = \frac{\text{transport time scale}} {\text{reaction time scale}} }. \]

Different chemical-engineering contexts use different specific forms, but the core purpose is the same: compare reaction speed with transport.

Diffusion Length

Over characteristic time \(t\), diffusion spreads over a scale of roughly

\[ \boxed{ L_D\sim\sqrt{Dt} }. \]

This scaling estimate helps determine whether a system can be treated as well mixed.

Well-Mixed Assumptions

Many ordinary differential-equation models assume uniform concentration throughout the reactor.

This approximation becomes questionable when reaction is faster than mixing, when strong gradients develop, or when geometry creates transport limitations.

Gas-Phase Models

For an ideal gas mixture,

\[ \boxed{ PV=nRT }. \]

Partial pressures satisfy

\[ \boxed{ P_i=y_iP } \]

where \(y_i\) is mole fraction under the ideal-mixture approximation.

Nonideal Systems

At high pressures, high concentrations, strong intermolecular interactions, or near phase transitions, ideal-gas and ideal-solution assumptions may fail.

More advanced models may use fugacity, activity coefficients, equations of state, or other corrections.

Acid-Base Models

Acid-base systems combine equilibrium relationships with conservation constraints.

For a simple acid dissociation

\[ HA\rightleftharpoons H^+ + A^-, \]

an idealized equilibrium expression is

\[ \boxed{ K_a= \frac{[H^+][A^-]}{[HA]} }. \]

Charge balance and total-material balance are also required.

Electrochemical Systems

Electrochemistry couples chemical species, charge transfer, electric potential, diffusion, and sometimes fluid flow.

Battery, fuel-cell, corrosion, and electrolysis models often require both reaction kinetics and transport equations.

Data and Parameter Estimation

Reaction parameters can be estimated from concentration-time data, conversion measurements, spectroscopy, calorimetry, pressure data, or other observations.

A least-squares objective may be

\[ \boxed{ J(\theta) = \sum_{i=1}^{n} \left[ y_i-\hat y(t_i;\theta) \right]^2 }. \]

Estimating a First-Order Rate Constant

For

\[ C_A(t)=C_{A,0}e^{-kt}, \]

taking logarithms gives

\[ \boxed{ \ln C_A =\ln C_{A,0}-kt }. \]

Under the model, the slope of \(\ln C_A\) versus \(t\) is \(-k\).

Model Selection

Several rate laws may fit a limited dataset. Model selection should consider more than numerical fit:

Identifiability

In multi-step mechanisms, different parameter combinations can sometimes produce nearly indistinguishable concentration curves.

A numerical optimizer may still report precise values, but the data may not actually contain enough information to identify each parameter separately.

Stiff Reaction Systems

Reaction networks can contain very fast and very slow processes at the same time. Such systems are often stiff.

Explicit time-stepping may then require extremely small time steps for stability even when the slow dynamics evolve over much longer periods.

Numerical Integration

For

\[ \frac{d\mathbf{C}}{dt} =\mathbf{f}(t,\mathbf{C}), \]

numerical solvers approximate the concentration trajectory through discrete time steps.

Solver choice should account for stiffness, accuracy, positivity, conservation, and event handling.

Positivity

Chemical concentrations should generally remain nonnegative. A numerical method that produces large negative concentrations has violated a physical constraint.

Small numerical undershoots may still require care, especially when those values feed nonlinear rate laws.

Numerical Conservation

A chemical solver should preserve elemental balances and other invariants to within justified numerical tolerance.

Drift in conserved quantities can reveal implementation mistakes or inadequate numerical methods.

Example: First-Order Decay

Suppose

\[ C_{A,0}=2.0\ \text{mol/L}, \qquad k=0.15\ \text{min}^{-1}. \]

After \(10\) minutes,

\[ C_A(10) =2.0e^{-0.15(10)} \] \[ =2.0e^{-1.5} \] \[ \boxed{ C_A(10)\approx0.446\ \text{mol/L} }. \]

Example: First-Order Half-Life

If

\[ k=0.20\ \text{h}^{-1}, \]

then

\[ t_{1/2} =\frac{\ln2}{0.20} \] \[ \boxed{ t_{1/2}\approx3.47\ \text{h} }. \]

Example: Arrhenius Ratio

If the same Arrhenius parameters apply at temperatures \(T_1\) and \(T_2\), then

\[ \boxed{ \ln\left(\frac{k_2}{k_1}\right) = -\frac{E_a}{R} \left( \frac{1}{T_2}-\frac{1}{T_1} \right) }. \]

This form allows comparison of rate constants without first determining the pre-exponential factor \(A\).

Example: CSTR Balance

Consider a constant-volume CSTR with first-order consumption

\[ R_A=-kC_A. \]

At steady state,

\[ 0=F C_{A,\text{in}}-FC_A-VkC_A. \]

Using \(\tau=V/F\),

\[ \boxed{ C_A= \frac{C_{A,\text{in}}}{1+k\tau} }. \]

Example: CSTR Conversion

Since

\[ X_A =1-\frac{C_A}{C_{A,\text{in}}}, \]

the first-order ideal CSTR result becomes

\[ \boxed{ X_A= \frac{k\tau}{1+k\tau} }. \]

This shows directly how conversion depends on the dimensionless group \(k\tau\).

Dimensionless Modeling

Suppose a first-order batch reaction satisfies

\[ \frac{dC}{dt}=-kC. \]

Define

\[ c=\frac{C}{C_0}, \qquad \tau=kt. \]

Then

\[ \boxed{ \frac{dc}{d\tau}=-c }. \]

The normalized dynamics no longer depend explicitly on \(k\) or \(C_0\).

Sensitivity Analysis

Sensitivity analysis examines how outputs respond to changes in parameters such as

Highly sensitive parameters deserve especially careful measurement and calibration.

Uncertainty

Chemical models may contain uncertainty in kinetic parameters, thermodynamic properties, transport coefficients, measurements, operating conditions, and model structure.

A model used for design or safety should propagate relevant uncertainty into outputs such as conversion, selectivity, temperature, pressure, or peak concentration.

Verification

Verification asks whether the equations have been implemented and solved correctly.

Useful checks include

Validation

Validation asks whether the model reproduces real chemical behavior well enough for the intended use.

Validation may compare measured and predicted concentration histories, conversion, selectivity, temperature, pressure, equilibrium composition, or spatial concentration profiles.

Common Error: Reading Reaction Order from Stoichiometry

The balanced overall reaction does not generally determine the kinetic rate law.

Reaction orders match stoichiometric coefficients automatically only under specific mechanistic assumptions, such as an elementary mass-action step.

Common Error: Ignoring Units of the Rate Constant

The units of \(k\) depend on the reaction order. For a first-order reaction,

\[ [k]=\text{time}^{-1}. \]

For other rate laws, the concentration dimensions must also be included so the reaction rate has correct units.

Common Error: Confusing Equilibrium with Completion

Equilibrium does not necessarily mean all reactants have been consumed. It means forward and reverse tendencies balance under the selected thermodynamic conditions.

Common Error: Treating Temperature as Fixed When It Is Not

In strongly exothermic or endothermic systems, reaction can change temperature enough to alter the rate constant materially.

A constant-temperature model may then miss the dominant feedback.

Common Error: Assuming Perfect Mixing

A well-mixed model can fail when reaction is faster than mixing or when diffusion, flow, geometry, or phase boundaries generate significant gradients.

In those cases, spatial transport must be modeled explicitly.

Common Error: Overfitting a Reaction Mechanism

A large mechanism with many adjustable rate constants may fit one dataset while remaining poorly identifiable and weakly predictive.

Mechanistic complexity should be supported by chemistry, data, and model purpose.

A Reliable Modeling Procedure

  1. Define the chemical question and the outputs that matter.
  2. List species, phases, reactions, and conserved quantities.
  3. Construct the stoichiometric relationships.
  4. Select kinetic and thermodynamic models appropriate to the mechanism.
  5. Add transport, flow, and energy balances when they materially affect behavior.
  6. Specify initial conditions, boundary conditions, and operating conditions.
  7. Estimate unknown parameters from suitable data.
  8. Verify units, conservation, positivity, and numerical convergence.
  9. Validate against measurements not used solely for calibration when possible.
  10. Perform sensitivity and uncertainty analysis before using the model for design or prediction.

Practice 1: Stoichiometric Change

For

\[ 2A+B\rightarrow3C, \]

the extent of reaction increases by \(0.40\) mol. Find the changes in \(n_A\), \(n_B\), and \(n_C\).

Practice 2: First-Order Concentration

A first-order reactant begins at \(1.50\,\text{mol/L}\) with \(k=0.10\,\text{min}^{-1}\). Find the concentration after \(8\) minutes.

Practice 3: First-Order Half-Life

Find the half-life when

\[ k=0.035\ \text{s}^{-1}. \]

Practice 4: Zero-Order Reaction

A zero-order reaction has

\[ C_{A,0}=0.80\ \text{mol/L}, \qquad k=0.050\ \text{mol/(L}\cdot\text{min)}. \]

Find \(C_A\) after \(6\) minutes.

Practice 5: Second-Order Reaction

For

\[ -\frac{dC_A}{dt}=kC_A^2, \]

let

\[ C_{A,0}=1.0\ \text{mol/L}, \qquad k=0.20\ \text{L/(mol}\cdot\text{min)}. \]

Find \(C_A\) after \(5\) minutes.

Practice 6: CSTR Residence Time

A reactor has volume \(V=2.0\,\text{m}^3\) and volumetric flow \(F=0.50\,\text{m}^3/\text{min}\). Find the residence time.

Practice 7: First-Order CSTR Conversion

For an ideal first-order CSTR,

\[ X_A=\frac{k\tau}{1+k\tau}. \]

Find conversion when \(k=0.25\,\text{min}^{-1}\) and \(\tau=4.0\,\text{min}\).

Practice 8: Michaelis-Menten Rate

Let

\[ V_{\max}=10\ \text{mmol/min}, \qquad K_M=2.0\ \text{mmol/L}, \qquad S=6.0\ \text{mmol/L}. \]

Find the reaction rate \(v\).

Practice 9: Diffusion Scale

Estimate the diffusion length after \(100\) seconds if

\[ D=1.0\times10^{-9}\ \text{m}^2/\text{s} \]

using \(L_D\sim\sqrt{Dt}\).

Practice 10: Choose the Better Model

A strongly exothermic reaction occurs in a poorly mixed vessel, and measured temperatures vary substantially across the reactor. Explain why a single, isothermal, well-mixed kinetic equation may be inadequate and identify the modeling effects that should be added.

Practice Answers

Practice 1

The signed stoichiometric coefficients are

\[ \nu_A=-2,\qquad \nu_B=-1,\qquad \nu_C=3. \]

Therefore,

\[ \Delta n_A=(-2)(0.40)=\boxed{-0.80\ \text{mol}}, \] \[ \Delta n_B=(-1)(0.40)=\boxed{-0.40\ \text{mol}}, \] \[ \Delta n_C=(3)(0.40)=\boxed{1.20\ \text{mol}}. \]

Practice 2

\[ C_A(t)=C_{A,0}e^{-kt} \] \[ =1.50e^{-(0.10)(8)} \] \[ =1.50e^{-0.8} \] \[ \boxed{ C_A\approx0.674\ \text{mol/L} }. \]

Practice 3

\[ t_{1/2} =\frac{\ln2}{0.035} \] \[ \boxed{ t_{1/2}\approx19.8\ \text{s} }. \]

Practice 4

\[ C_A=C_{A,0}-kt \] \[ =0.80-(0.050)(6) \] \[ \boxed{ C_A=0.50\ \text{mol/L} }. \]

Practice 5

\[ \frac{1}{C_A} =\frac{1}{C_{A,0}}+kt \] \[ =1+(0.20)(5) =2. \]

Therefore,

\[ \boxed{ C_A=0.50\ \text{mol/L} }. \]

Practice 6

\[ \tau=\frac{V}{F} =\frac{2.0}{0.50} =\boxed{4.0\ \text{min}}. \]

Practice 7

\[ k\tau=(0.25)(4.0)=1. \] \[ X_A=\frac{1}{1+1} =\boxed{0.50}. \]

The predicted conversion is \(50\%\).

Practice 8

\[ v= \frac{V_{\max}S}{K_M+S} \] \[ =\frac{(10)(6)}{2+6} \] \[ \boxed{ v=7.5\ \text{mmol/min} }. \]

Practice 9

\[ L_D\sim\sqrt{Dt} \] \[ =\sqrt{ (1.0\times10^{-9})(100) } \] \[ =\sqrt{1.0\times10^{-7}} \] \[ \boxed{ L_D\approx3.16\times10^{-4}\ \text{m} } \]

or about \(0.316\,\text{mm}\).

Practice 10

The assumptions of uniform concentration and constant temperature are not supported by the observed system. The model should include mixing or spatial transport, such as advection and diffusion, together with an energy balance that couples temperature to reaction rate through temperature-dependent kinetics. Depending on geometry, this may require a spatial reaction-transport model rather than a single ordinary differential equation.

What to Remember

Where We Go Next

Chemical reaction models show how conservation laws, nonlinear rates, equilibrium, transport, and feedback combine to create system behavior. The same mathematical tools extend naturally into economic and operational systems, where stocks, flows, prices, inventories, demand, risk, and decision rules interact over time.

The next lesson, Economics, Finance, and Operations, develops applied models of growth, markets, cash flows, inventories, queues, optimization, risk, and operational decision systems.