Sources and References

• Senturia, Microsystem Design, Springer.
• Langtangen and Linge, Finite Difference Computing with PDEs: A Modern Software Approach, Springer.
• Reddy, An Introduction to the Finite Element Method, 4th ed., McGraw-Hill.
• Quarteroni, Sacco, and Saleri, Numerical Mathematics, 2nd ed., Springer.
• LeVeque, Finite Volume Methods for Hyperbolic Problems, Cambridge University Press.

Chapter 1: Modelling and Simulation for Nanosystems

Computer-aided design of micro and nanosystems combines physical modelling, numerical analysis, and software tools into a design loop: formulate the physics, discretize, solve, interpret, iterate. The engineer’s goal is not a beautiful simulation but an efficient design; simulation is a means to that end.

1.1 Why Simulate

Nanosystems are expensive to fabricate and to characterize. Simulation predicts performance before commitment of a mask set or a process run, identifies sensitivities for robust design, and interprets characterization data. It is indispensable but not infallible: models are only as good as their physics and their parameters.

1.2 Fidelity and Cost

Simulation fidelity scales with cost. A closed-form analytical model may be instant but limited; a reduced-order model captures dominant dynamics quickly; a three-dimensional finite-element analysis handles geometric complexity at greater cost; a coupled-multiphysics simulation with nonlinearities and many coupled variables is the most expensive. The engineer selects the coarsest adequate model.

Chapter 2: Lumped Versus Distributed Modelling

2.1 Lumped Models

A lumped model represents distributed physics by a finite number of nodes connected by idealized elements. Electrical circuits, spring–mass–damper mechanical systems, and thermal RC networks are lumped. The approximation is accurate when characteristic lengths of interest are larger than the domain of interest (spatially uniform variables) or when the system has been deliberately designed to operate in a regime where distributed effects are negligible.

A lumped MEMS cantilever with stiffness \( k \), mass \( m \), and damping \( c \) obeys

and couples to an electrostatic actuator with force

2.2 Distributed Models

When variables vary substantially over the system, partial differential equations replace ODEs. Beam bending, heat conduction, and diffusion all take PDE form. The extra generality costs computational effort and demands discretization.

2.3 Reduced-Order Models

Distributed models can be projected onto a small basis — Galerkin with dominant mode shapes, proper orthogonal decomposition of simulation snapshots — producing reduced-order models that preserve the key dynamics at a fraction of the cost. Reduced-order models serve for system-level simulation in Simulink or Verilog-A environments and for real-time control-design iteration.

Chapter 3: Differential Equations and Constitutive Relations

3.1 The Standard PDE Families

Elliptic equations — Laplace, Poisson — describe equilibrium: electrostatics \( \nabla^2 \phi = -\rho/\varepsilon \), steady-state temperature, steady diffusion. Parabolic equations — heat, diffusion — describe transient relaxation. Hyperbolic equations — wave, advection — describe propagation.

3.2 Coupled Problems in Microsystems

Most microsystems are multi-physics. An electrostatically actuated MEMS switch couples electrostatics (distribution of electric field), elasticity (displacement of the movable electrode), and sometimes squeeze-film gas damping. A micro-fluidic sensor couples fluid flow (Navier–Stokes), species transport (advection–diffusion), and surface binding kinetics. A piezoelectric resonator couples elasticity and electrical response through the piezoelectric tensor.

3.3 Constitutive Relations

Constitutive relations close the PDEs by linking fluxes to state variables. Linear elasticity relates stress to strain; Ohm’s law relates current to field; Fick’s law relates diffusive flux to concentration gradient; Fourier’s law relates heat flux to temperature gradient. At the nanoscale, classical constitutive relations may need modification: size-dependent elasticity, non-Fourier heat conduction, slip flow in confined gases.

Chapter 4: Numerical Solution Methods

4.1 Finite Difference

Finite-difference methods discretize derivatives on a regular grid. Forward, backward, and central differences give first-order or second-order accuracy. Explicit time-stepping is simple but restricted by CFL or diffusion stability limits; implicit schemes are unconditionally stable but require solving linear systems each step. Finite-difference codes are easy to write and adequate for many rectangular-geometry problems in MEMS.

4.2 Finite Element

Finite-element methods discretize the weak form of a PDE on an unstructured mesh. Shape functions (linear, quadratic, higher-order) interpolate within each element; assembly produces a global linear (or, for nonlinear problems, a sequence of linearized) system. FEM handles arbitrary geometry, material heterogeneity, and coupled fields; it is the default in commercial multiphysics tools.

4.3 Finite Volume

Finite-volume methods, dominant in CFD, enforce conservation on discrete control volumes. Flux interpolation between cells determines accuracy and stability. They extend naturally to nanofluidic problems, including multiphase flow and reactive transport.

4.4 Boundary Element and Meshless Methods

Boundary-element methods discretize only the boundary of the domain — excellent for exterior problems and for near-field electromagnetics. Meshless methods avoid the meshing step at the cost of more expensive assembly; they are used for large-deformation solid mechanics and nanoscale simulations where topology changes.

Chapter 5: Solvers for Nonlinear and Coupled Problems

5.1 Newton’s Method

Nonlinear algebraic equations \( \mathbf{F}(\mathbf{x}) = 0 \) are solved by Newton iteration,

with Jacobian \( J = \partial \mathbf{F}/\partial \mathbf{x} \). Line search, trust-region modifications, and quasi-Newton approximations extend robustness and efficiency. Convergence near a solution is quadratic; global convergence depends on the initial guess.

5.2 Coupling Strategies

Coupled multiphysics solvers use either monolithic or segregated strategies. Monolithic approaches assemble all equations into a single system and solve together; they are robust to strong coupling but may be expensive. Segregated approaches alternate between sub-problems, applying Picard or Newton iterations; they scale better for loosely coupled systems but may fail for strong nonlinear coupling such as pull-in instabilities.

5.3 Time Integration

Time-dependent problems need integration schemes suited to the stiffness and accuracy demands. Explicit Runge–Kutta for mild stiffness; implicit BDF or Crank–Nicolson for stiff problems. For oscillatory problems of interest in MEMS resonators, symplectic integrators preserve energy to acceptable drift over long simulations.

Chapter 6: Boundary Conditions and Modelling of Nanoscale Physics

6.1 Boundary Conditions

Dirichlet (prescribed value), Neumann (prescribed flux), and Robin (mixed) conditions cover the typical cases. Symmetry and periodic conditions exploit geometric regularity. Special conditions — contact, radiation, absorbing boundaries for wave problems — require additional care.

6.2 Molecular-Continuum Coupling

At some nanoscale interfaces the continuum description breaks down: a fluid in a channel of a few molecular diameters is no longer well described by Navier–Stokes with no-slip walls; heat conduction in a thin film with phonon mean free path comparable to thickness is no longer Fourier. Hybrid molecular–continuum schemes couple atomistic regions near interfaces with continuum domains in the bulk, enabling simulation across scales.

6.3 Fluidic Systems

Microfluidic and nanofluidic modelling solves the Navier–Stokes equations in complex channel geometries, often coupled to species transport and electrokinetic forcing. Electrokinetic flow under applied electric fields drives separations and pumps without moving parts. Slip boundary conditions become important as channel dimensions approach the molecular mean free path.

6.4 Optical, Thermal, Mechanical Micro- and Nanostructures

Optical MEMS require electromagnetic simulation — FDTD, frequency-domain FEM, or beam propagation — coupled to mechanical deformation and sometimes to thermal tuning. Thermal MEMS demand conjugate heat transfer including radiation. Mechanical nanostructures benefit from size-dependent elasticity where classical elasticity overpredicts stiffness at small scales.

Chapter 7: From Model to Design

7.1 Verification and Validation

Verification checks that the solver solves the equations correctly (mesh convergence, method-of-manufactured-solutions tests). Validation checks the equations against experiment. A design supported by verified, validated simulation is defensible; one supported by an unvalidated simulation is a bet.

7.2 Parameter Studies and Optimization

Simulation is most valuable for exploring design spaces. Parameter sweeps reveal sensitivities; design-of-experiments planners (factorial, Latin-hypercube) efficiently sample multidimensional spaces; optimization drives designs to best combinations of objectives. Surrogate models (polynomial regression, Gaussian processes, neural networks) accelerate expensive multiphysics loops.

7.3 Workflow and Reproducibility

A reproducible simulation workflow documents geometry, mesh, material parameters, boundary conditions, solver settings, and post-processing choices. Scripts that automate the workflow eliminate transcription error, enable version control, and let collaborators repeat the analysis. In an academic course and in industry alike, these habits separate trustworthy results from unrepeatable ones.